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摘 要 : 以 塔里木 盆地 南 缘 的 策 勤 绿 洲 为 研究 区 , 基 


据 , 分 析 了 1990 一 2018 年 策 勒 绿洲 土地 利用 / 履 被 变化 特 和 


于 4 期 土地 利用 和 遥感 影像 ,以 及 气象 .水文 和 社会 经 济 数 
F。 结 果 表 明 :(1) 策 勒 绿洲 各 土地 利用 类 型 发 生 明 显 转 


变 , 农 田 面 积 向 四 周 扩 增 ,林地 分 布 区 域 发 生 大 幅度 转移 ,不同 覆盖 程度 的 草地 之 间 转 化 频繁 ,(2) 绿洲 共有 138.41 
km 的 面积 发 生 转 变 , 占 总 面积 的 53.85% ,各 类 型 的 单一 动态 度 依 次 为 :高 覆盖 草地 > 农田 > 低 覆 盖 草 地 > 建设 用 地 > 
中 覆盖 草地 > 未 利用 土地 > 林地 ,双向 动态 度 依 次 为 :高 覆盖 草地 > 农田 > 低 履 盖 草 地 > 建设 用 地 > 林地 > 中 覆盖 草地 > 


A 
人 口 降水 以 及 径流 量 。 


c— 


] 土 地 。(3) 绿洲 林 草 植被 指数 低 于 0, 表 明 林 地 和 草地 


直 处 于 退化 的 状态 ;绿洲 土地 利用 变化 驱动 力主 要 为 
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土地 利用 / 窗 盖 变化 是 分 析 经 济 .社会 和 上 自然 环 
境 交 织 作用 的 核心 ,在 全 球 环境 变化 过 程 中 扮演 着 
至 关 重 要 的 角色 。 土 地 利用 是 特定 人 地 关系 下 
多 种 因素 共同 作用 的 结果 ,土地 作为 陆地 生态 系统 
的 基础 , 即 是 人 类 生存 的 依托 ,又 是 人 类 与 环境 相 
互 作 用 的 主要 表现 形式 ”。 随 着 现代 社会 的 进步 和 
生产 力 的 大 力 发 展 ,人 们 对 土地 的 利用 方式 不 断 转 
变 、 利 用 强度 日 益 增 加 ,导致 土地 利用 类 型 发 生 巨 
大 变化 “"。 目 前 ,国内 外 学 者 围绕 土地 利用 / 覆 被 变 
化 的 不 同时 间 和 空间 尺度 ,借助 遥感 和 地 理 信息 系 
统 技术 平台 ,结合 多 种 理论 和 方法 ,运用 数学 模型 、 
结构 模型 等 ,针对 土地 利用 的 时 空 过 程 .模拟 预测 、 
动力 学 机 制 和 景观 格局 展开 广泛 的 研究 一 。Li 等 ” 
学 者 对 土地 利用 模拟 技术 的 空间 分 辨 率 作 了 改进 ， 
为 研究 全 球 土地 利用 变化 提供 了 一 种 更 为 准确 的 
方法 ; 刘 纪 元 等 ”提出 的 通过 遥感 图 像 和 人 机 交互 
解 译 法 为 分 析 土 地 利用 变化 提出 了 一 种 新 思路 。 
不 少 学 者 以 生境 质量 .景观 结构 .土地 利用 强度 、 风 
险 评估 等 为 切入 点 对 土地 利用 变化 的 特征 进行 了 
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详细 的 描述 ””。 大 量 研究 表明 ,区域 变化 综合 效应 
可 能 会 带 来 对 流域 尺度 乃至 全 球 尺度 的 影响 "， 
此 研究 区 域 尺 度 LUCC 时 空 特征 及 驱动 力 是 研究 大 
尺度 变化 力 至 全 球 变 化 演变 过 程 的 基础 。 

在 区 域 尺 度 方面 以 绿洲 土地 利用 变化 的 研究 
最 为 典型 ,中 国 的 绿洲 基本 都 分 布 在 西北 地 区 ,如 
新 疆 、 甘 肃 、 宁 夏 和 内 蒙古 等 地 区 ,其 中 大 部 分 绿洲 
主要 分 布 在 新 疆 , 其 绿洲 承载 了 新 疆 90% 以 上 的 耕 
地 、 人 口 及 经 济 生产 总 值 。 然 而 ,绿洲 生态 环境 极 
其 脆弱 ,人 类 活动 的 微弱 的 变化 就 会 造成 绿洲 生境 
破碎 的 风险 增加 ,进而 影响 绿洲 的 区 域 生态 及 
水 文 过 程 。 目 前 , 受 人 类 活动 影响 ,干旱 区 绿洲 的 
土地 利用 / 窗 盖 格局 正在 发 生 快速 变化 ,并 且 这 种 变 
化 对 绿洲 脆弱 生态 环境 平衡 产生 巨大 影响 ”。 为 
改善 不 合理 的 土地 利用 方式 ,提高 土地 利用 可 持续 
发 展 ,学 者 们 从 数量 变化 .空间 变化 .土地 利用 集约 
程度 变化 .景观 格局 及 动态 度 方面 分 析 绿 洲 土 地 利 
用 变化 的 特征 ,并 应 用 相关 分 析 和 主 成 分 分 析 对 影 
响 绿洲 土地 利用 变化 的 驱动 因素 进行 探讨 ,以 改善 
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绿洲 的 生态 环境 ,提出 绿洲 可 持续 性 发 展 的 相应 对 
gn 

近 几 十 年 来 , 策 勒 绿洲 在 人 口 增长 和 经 济 发 展 
的 驱动 下 ,不 断 开采 抽水 井 . 开 以 天 然 植 被 为 驮 地 ， 
使 得 绿洲 耕地 面积 及 灌溉 量 大 幅 增加 ,绿洲 景观 格 
局 不 断 发 生变 化 ,绿洲 退化 风险 程度 加 剧 。 因 此 ， 
本 文 以 新 疆 塔里木 盆地 南 缘 的 策 勒 绿洲 为 研究 区 ， 
通过 分 析 长 时 间 尺 度 绿洲 地 表 覆 被 时 空 变化 及 其 
了 驱动 机 制 , 为 保护 绿洲 生态 环境 ,制定 科学 土地 空 
间 管 理 策略 ,实施 可 持续 发 展 目标 提供 重要 的 科学 
依据 。 


1 研究 区 概况 与 研究 方法 


1.1 研究 区 概况 

策 勒 绿洲 地 处 塔里木 盆地 南 缘 (图 1) ,隶属 于 新 
疆 和 田地 区 策 勒 县 ,其 边界 以 绿洲 分 布 范 围 进行 界 
定 ,地理 坐标 为 35°18’~39°18’N,80°03’~82°10’E。 
其 地 势 呈 南 高 北 低 ,海拔 在 1296.5~1370.5 m 之 间 ， 
中 部 发 育 山药 倾斜 平原 。 流 经 绿洲 的 策 勒 河 年 均 
径流 量 为 1.28xl105m。 绿 洲 地 区 气候 为 暖 温带 干旱 
荒漠 型 ,气温 年 均值 为 11.9 % ,多 年 极端 最 高 气温 
为 42.0% ,极端 最 低 气温 为 -23.9 C ,年 平均 降水 量 
仅 为 35.1 mm, 蒸 发 量 接近 2595.3 mm, Z&Z Et ir e 
超过 降水 量 ,全 年 无 霜 期 在 230 d 左 右 , 全 年 日 照 时 
数 在 2686h 以 上 。 


À 


um 农田 
umm 林地 

mm 水 体 

= 草地 

c 建设 用 地 
E 未 利用 土地 


图 1 策 勒 绿洲 位 置 示意 图 


Fig. 1 Overview of Cele oasis 
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1.2 数据 来 源 

采用 来 自 中 国 科学 院 资 源 环 境 数 据 云 平台 
(http://www.resde.cn/) 提 供 的 1990 年 、2000 年 .2010 
年 .2018 年 的 策 勒 绿洲 土地 利用 数据 ,其 分 辨 率 均 为 
30 mx30 m, 土 地 利用 划分 为 耕地 、 林 地 水域. 建设 
用 地 草地、 未 利用 土地 6 个 一 级 土地 利用 类 型 , 草 
地 进一步 划分 为 高 履 盖 草地 中 和 窗 盖 草地 和 低 履 盖 
草地 3 个 二 级 土地 利用 类 型 ,该 数据 通过 了 质量 检查 
和 精度 检验 ,土地 利用 一 级 类 型 综合 评价 精度 达到 
94.3% 以 上 ,二 级 类 型 分 类 综合 精度 达 91.2% 以 上 2 ， 
且 该 数据 目前 已 被 广泛 应 用 于 全 国 范围 内 不 同 空间 
时 间 尺 度 的 土地 利用 /土地 覆盖 变化 监测 9。 

还 感 影像 来 自 美国 地 质 勘 探 局 (http://earthex- 
plorer.usgs.gov/) 用 于 计算 植被 的 NDVI, 其 轨道 号 为 
146/34 ,选取 当日 天 气 晴朗 , 云 量 较 低 ,像素 分 辩 率 
均 为 30 mx30 m, 采 用 与 土地 利用 数据 年 份 相应 的 4 
期 遥感 影像 ,为 反映 绿洲 植被 真实 的 生长 状况 ,下 
载 了 生长 季 期 间 的 遥感 影像 ,分别 为 1990 年 Land- 
sat5 TM(7 月 18 号 ) 2000 4E Landsat7 ETM+(8 月 5 
号 )、2010 年 Landsat5 TM(7 月 8 号 )、2018 4E Land- 
sat8 OLI(7 月 14 号 ) 共 4 期 卫星 影像 。 对 4 期 遥感 影 
像 首先 运用 ENVI 软 件 进 行 辐射 定 标 .大气 校正 .图 
像 裁剪 等 预 处 理 ,再 根据 不 同 遥 感 影像 的 传感器 选 
择 相 应 波段 计算 NDVI 并 去 除 异 常 值 ,最 终 得 到 4 期 
策 勒 绿洲 NDVI 空间 分 布 数据 。 以 这 4 期 愧 感 影像 
为 基准 ,利用 ArcGIS 在 研究 区 范围 内 随机 生成 100 
个 点 ,这 100 个 随机 点 覆盖 了 研究 区 的 所 有 土地 利 
用 类 型 ,能够 充分 的 代表 各 土地 利用 类 型 ,对 来 自 
中 国 科 学 院 资源 环境 中 心 的 1990 年 .2000 年 .2010 
年 和 2018 年 土地 利用 数据 进行 验证 ,在 验证 过 程 中 
对 于 不 确定 土地 利用 类 型 结合 同期 谷歌 历史 高 清 
影像 进行 地 物 识 别 ,验证 结果 显示 ,4 期 土地 利用 数 
据 一 级 分 类 精度 分 别 为 94% .95% 96% 95%. 

本 文 用 到 的 年 尺度 气温 、 降 水 和 蒸发 量 数据 来 
自 于 策 勒 气象 站 ,径流 量 数据 来 自 策 勒 水 文 站 ,这 
些 数据 之 前 有 学 者 进行 过 分 析 研 究 ””, 在 本 文中 
主要 是 用 于 分 析 较 长 时 间 尺 度 上 对 策 勒 土地 利用 
的 影响 。 
1.3 研究 方法 

首先 基于 策 勒 绿洲 4 期 土地 利用 数据 分 析 该 绿 
洲 各 土地 利用 类 型 的 时 空 分 布 特征 ,再 通过 土地 利 
用 转移 矩阵 计算 各 土地 利用 类 型 数量 上 的 变化 ,其 
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次 根据 林 草 植被 指数 判断 转化 过 程 中 绿洲 林地 和 
草地 的 增 减 状况 ,以 及 利用 动态 度 分 析 各 土地 利用 
类 型 的 变化 程度 ,最 后 通过 MK 趋势 分 析 、 主 成 分 分 
析 及 多 元 线性 回归 模型 等 方法 ,分 析 策 勒 绿洲 近 几 
十 年 来 土地 利用 类 型 变化 的 驱动 因素 。 

1.3.1 土地 利用 转移 矩阵 采用 马尔 可 夫 转 移 矩 阵 
模型 对 策 勒 绿洲 各 土地 利用 类 型 之 间 数 量 上 的 转 
化 情况 进行 描述 ,其 数学 表达 式 ” 如 下 所 示 : 


Pu Pa $9 Pin 
Ps Pa Py + P, (1) 
m Ps Ut B mn 


RP: n 为 土地 利用 类 型 ; P, 为 土地 利用 类 型 m 
转化 为 土地 利用 类 型 n 的 转移 概率 , P,, 须 同时 满 
足 2 个 条 件 :0< P, «1, M P, 71. 
1.3.2 林 草 植被 变化 指数 ”该 指数 反映 了 某 一 个 地 
区 植被 覆盖 变化 程度 的 重要 指标 ,文中 用 于 分 析 绿 
洲 林地 和 草地 的 增 减 情况 ,计算 公式 如 下 ™. 

T=F, -PF, (2) 
式 中 :了 是 林 草 植被 变化 指数 , 取 值 -1~1; Fan Fo 
4] ile tl , 刀 年 林 草 用 地 面积 比重 。 若 工 值 为 正 ， 
表明 区 域 林 草 植被 增加 , 且 了 值 越 大 ,说 明 增加 的 速 
度 越 快 ; 若 了 值 为 负 , 则 反之 。 
1.3.3 土地 利用 类 型 动态 度 单一 土地 利用 类 型 动 
态度 表达 了 一 定时 间 范 围 内 各 土地 利用 类 型 的 数 
量变 化 情况 ,其 值 反 映 绿洲 某 类 土地 利用 的 变化 幅 
度 和 速度 ,其 表达 式 为 : 


U,-U 
Ke e x ex 100% 


(3) 


AP: K 为 研究 时 段 内 某 一 土地 利用 类 型 动态 度 ; 
U,、U, 分 别 为 研究 期 初 和 研究 期 末 某 一 种 土地 利 
用 类 型 的 数量 ; 7 为 研究 时 段 长 。 符 K 值 为 正 数 ， 
则 土地 利用 类 型 的 面积 倾向 于 增加 的 趋势 , 千 K 为 
负 值 , 则 反之 ,的 绝对 值 大 小 表示 某 种 土地 利用 
类 型 转变 的 快慢 。 

双向 土地 利用 度 态 度 表 示 某 一 土地 利用 类 型 
扩张 与 退缩 双向 动态 变化 率 , 可 以 反映 各 土地 利用 
类 型 之 间 相 互 转化 的 剧烈 程度 ,其 表达 式 为 ”: 

(A U,+AU)) 
U, 
式 中 :天 为 第 ;类 土地 利用 类 型 的 双向 变化 率 ; AU, 
表示 由 第 i 类 土地 利用 类 型 转移 到 其 他 土地 利用 类 


K= x $ x 100% (4) 


型 的 面积 ; A UV 为 同时 期 其 他 土地 利用 类 型 转变 为 
第 i 类 土地 利用 类 型 的 面积 ; U, 为 研究 初期 第 i 类 
土地 利用 类 型 的 面积 。 
1.3.4 计算 归 一 化 植被 指数 (NDVI) 植被 指数 是 
通过 卫星 某 些 波段 对 植被 敏感 的 特性 ,将 这 些 光谱 
数据 结合 到 一 起 所 形成 的 图 像 植 被 指数 ,被 用 来 分 
析 绿 洲 植被 生长 状况 ,其 计算 公式 如 下 ”. 

NDVI =(NIR - R)/(NIR + R) (5) 
式 中 : R 是 可 见 光 红 波 段 ; NIR 是 近 红 外 波段 ;其 中 
TM 传感器 的 红 光 波段 为 B3; 近 红外 波段 为 B4; 
ETM+ 与 0LT 传 感 器 的 红 光 波段 为 B4, 近 红外 波段 
为 B5; 计 算 结果 将 得 到 一 幅 灰 度 值 表 示 的 单 波段 影 
像 ,其 像 元 值 范 围 一 般 在 -1< NDVI& 1 。 
1.3.5 M-K 趋 势 分 析 Mann-kendall 检验 对 样本 的 
分 布 特殊 要 求 ,数据 独立 即 可 ,基本 不 受 离 群 值 的 
影响 ,文中 用 于 分 析 各 驱动 因素 的 趋势 变化 特 
征 ,在 检验 中 标准 的 正 态 统计 变量 z 计算 公式 如 
Fr: 


S-1 


x a s»0 
Var(S) 
z= 0 if S=0 (6) 
S*1  ifg«0 
Var(S) 


S 为 正 态 分 布 ,其 均值 为 0, 检 验方 法 为 在 给 定 
的 o 置信 水 平 上 ,如 果 |z|> Zal > 则 原 假 设 是 不 可 
接受 的 , 即 在 a 置信 水 平 上 ,时 间 序 列 数据 存在 明 
显 的 上 升 或 者 下 降 趋势 。 对 于 统计 变量 z>0 时 ,是 
向 上 趋势 ,反之 , 则 是 下 降 趋 势 。 
13.6 多 元 线性 回归 ”表示 多 个 自 变量 与 因 变 量 之 
间 是 线性 关系 ,本 文中 用 于 拟 合 各 驱动 因素 与 土地 
利用 类 型 变化 的 关系 ,本 文中 将 土地 利用 类 型 设 为 
因 变 量 7, 将 各 驱动 因素 设 为 自 变 量 X ,计算 公式 
如 下 : 


Y-ayjta,X;*ta,X,* +a,X, (7) 
式 中 E 了 表示 因 变 量 ; A As UV, 表示 自 变 量 ; 
ak 表示 回归 系数 。 


2 结果 与 分 析 


2.1 绿洲 土地 利用 时 空 变化 
2.1.1 土地 利用 分 布 与 变化 由 图 2 所 示 , 策 勒 绿 
洲 4 期 土地 利用 分 布 表 现 为 农田 分 布 在 绿洲 内 部 ， 


d, d», .. 
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(a) 1990 年 


um 农田 
um 林地 
um 水 体 
= 建设 用 地 
= 高 覆盖 草地 
LIEU OL 
um 低 覆 盖 草 地 
= 未 利用 土地 LJ 


>z 


(c) 20104 


um 农田 
um 林地 
m 水 体 
= 建设 用 地 
= 高 覆盖 草地 
ma 中 覆盖 草地 
um 低 覆 盖 草 地 
= 未 利用 土地 
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(b) 2000 年 


>z 


um 农田 
um 林地 
um 水 体 
= 建设 用 地 
cc 高 覆盖 草地 
mm 中 覆盖 草地 
um 低 覆 盖 草 地 
e 未 利用 土地 NET 


(d) 2018 年 


- 


m 未 利用 土地 


图 2 策 勒 绿洲 4 期 土地 利用 分 布 情况 


Fig.2 Distribution of land use in Cele oasis in four periods 


占据 绿洲 的 面积 比重 越 来 越 大 ;不 同 覆 盖 程 度 的 草 
地 均 分 布 在 农田 的 外 围 ; 林 地 在 1990 41 2000 年 
4j i XE 4& EH PIER] , fk 2010 4E fI 2018 年 绿洲 零散 的 
分 布 在 农田 内 部 ,未 利用 土地 位 于 绿洲 边缘 ,面积 


料 的 需求 而 进行 的 伐木 砍 枝 ; (2) 策 勒 绿洲 东北 部 
建立 较为 完善 的 渠 系 结构 , 随 着 农田 在 该 方向 扩 增 
林地 面积 减少 ;(3) 策 勒 河上 游 农田 灌溉 量 增加 , 河 
水 无 法 到 达 下 游 使 植被 生长 受 限 以 及 绿洲 边缘 地 


不 断 的 缩减 ,建设 用 地 分 布 在 绿洲 内 部 ,整体 而 言 ， 
策 勒 绿洲 各 土地 利用 类 型 的 分 布 位 置 呈现 出 空间 
集约 性 ,并 且 各 土地 利用 类 型 的 分 布 从 内 向 外 扩 增 
延伸 。 

1990—2018 年 间 策 勒 绿洲 各 土地 利用 类 型 的 
面积 和 空间 分 布 发 生 显著 变化 (图 3) ,农田 分 布 面 
积 明 显 增 大 ,逐渐 成 为 占据 绿洲 的 主要 土地 利用 类 
型 ,增加 区 域 集中 在 外 围 , 尤 其 在 东北 部 明显 扩 
增 。 林 地 的 分 布 位 置 发 生 大 转变 ,1990 年 主要 分 布 
在 绿洲 东北 部 ,2000 年 集中 分 布 在 西北 部 ,2010 年 
和 2018 年 林地 均 交 错 分 布 在 农田 内 部 ,最 终 表现 为 
林地 面积 大 幅度 减少 ,新 增 部 分 较 少 且 分 布 零散 。 
不 同 覆 盖 程 度 的 草地 分 布 范 围 逐 渐 向 绿洲 边缘 缩 
减 ,上 且 不 同 草地 类 型 内 部 相互 发 生 转 变 。 究 其 林 草 
地 退化 的 原因 :(1) 居民 对 燃料 . 居 房 建材 和 家 冀 饲 


区 整体 地 下 水 位 下 降 ,进而 植被 草地 生长 衰退 以 至 
消 记 5。 建设 用 地 面积 在 2000 年 明显 增加 ,在 2010 
年 表现 为 中 心 城镇 的 扩 增 ,在 2018 年 表现 为 在 绿洲 
东南 部 新 增 , 但 增加 幅度 较 小 。 绿 洲 未 利用 土地 减 
少 部 分 明显 ,新 增 较 少 。 

2.1.2 绿洲 土地 利用 转移 情况 ”利用 ArcGIS 空 间 镭 
加 工具 得 到 各 种 土地 利用 类 型 之 间 上 有 具体 的 转化 情 
况 ,一 方面 可 以 分 析 不 同 土地 利用 类 型 空间 上 的 变 
化 过 程 ; 另 一 方面 可 以 分 析 不 同 土地 利用 类 型 数量 
上 的 稳定 性 特征 ,计算 结果 如 表 1 所 示 ,研究 时 段 内 
策 勒 绿洲 各 土地 利用 类 型 变化 表现 为 :农田 .居民 
建设 用 地 面积 增加 ,林地 、 未 利用 土地 、 高 覆盖 草 
地 .中 履 盖 草地 和 低 覆 盖 草 地 面积 减少 。 最 终 绿洲 
共有 138.41 km 的 面积 发 生 转 变 , 占 总 面积 的 
53.8596 ,其 中 耕地 为 主要 的 转 和 者 、 林 草 地 .未 利用 
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(a) 农田 À (b) 林地 À 


c 减少 部 分 
um 未 变 部 分 
um 增加 部 分 


(c) 高 覆盖 草地 


>z 


m AUD 
== 增加 部 分 0 2km 


(e) 低 覆 盖 草 地 


>z 


= 减少 部 分 
um 未 变 部 分 
mm 增加 部 分 


= 减少 部 分 
um 未 变 部 分 
um 增加 部 分 


mn FÉ 
EE 增加 部 分 0  3km 


(d) 中 覆盖 草地 


[zz] FÉ 
um 增加 部 分 o gim 


(建设 用 地 


m 


= 减少 部 分 
um 未 变 部 分 
um 增加 部 分 


图 3 策 勒 绿洲 土地 利用 类 型 的 空间 变化 


Fig.3 Spatial variation of land use type of Cele oasis 
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表 1 1990—2018 FRPR JU T- HR FH SETS AB EE 
Tab.1 Land use torque matrix of Cele oasis from 1990 to 2018 /km 

土地 利用 类 型 农田 林地 高 覆盖 植被 。 中 覆盖 植被 — 低 覆 盖 植 被 。 建设 用 地 ”未 利用 土地 总 计 

KI 56.55 5.02 0.24 0.02 1.68 0.10 63.61 
林地 33.56 1.24 7.05 1.45 0.02 0.13 0.05 43.50 
高 覆盖 植被 4.20 2.32 0.62 7.15 

中 覆盖 植被 12.75 16.04 17.79 1.38 2.55 50.51 

低 覆 盖 植 被 1.58 0.00 1.82 0.21 0.13 3.74 
建设 用 地 2.62 0.37 2.41 5.40 
未 利用 土地 23.15 0.22 2.36 11.19 4.58 3.52 38.09 83.11 
总 计 134.41 6.84 27.78 33.11 6.21 7.75 40.92 257.02 
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土地 为 主要 的 转 出 者 。 具 体 表现 为 :林地 77% 转 化 
为 农田 沙漠 28% 转 化 为 农田 高 覆盖 草 58% 转 化 为 
农田 中 覆盖 草地 25% 转 化 为 农田 . 低 覆 盖 草 地 42% 
转化 为 农田 以 及 未 利用 土地 27.85% 转 化 为 农田 , 农 
田 面积 扩 增 了 70.80 km ,成 为 策 勒 绿洲 面积 最 大 的 
土地 利用 类 型 。 建 设 用 地 的 增加 主要 来 自 未 利用 
土地 , 占 比 为 45.44%。 绿 洲 林地 和 草地 增长 面积 > 
52.37 km ,退化 的 面积 为 83.34 km ,最 终 退 化 30.97 
km 。 再 者 通过 计算 得 到 1990 一 2018 年 的 该 绿洲 林 
草 植被 变化 指数 均 小 于 0, 其 中 1990 一 2000 年 该 指 
数 为 -2.23% , 2000—2010 年 该 指数 为 -1.75% , 
2010—2018 年 该 指数 为 -8.07% ,在 2010—2018 年 
间 的 林 草 植被 指数 最 小 ,表明 在 这 个 时 间 段 策 勒 绿 
洲 的 林 草 植被 减少 的 速度 加 快 ,2000 一 2010 年 间 的 
林 草 植被 指数 最 大 ,表明 在 这 个 时 间 段 策 勒 绿洲 林 
草地 减少 ,但 减少 速度 在 整个 研究 时 间 段 最 小 。 在 
整个 研究 时 间 段 林 草 指数 均 为 -12.05% ,说 明 在 
1990 一 2018 年 , 策 勒 绿洲 植被 (林地 和 草地 ) 一 直 处 
于 减少 的 状态 。 


2.2 绿洲 土地 利用 动态 度 变 化 分 析 

为 进一步 分 析 各 土地 利用 类 型 转化 的 速率 , 计 
算 了 绿洲 单一 动态 度 和 双向 动态 度 如 下 表 2, 在 
1990—2018 年 间 各 土地 利用 类 型 的 动态 度 变 化 表 
现 为 :农田 高 覆盖 植被 低 覆 盖 植 被 和 建设 用 地 的 
单一 动态 度 和 双向 动态 度 均 大 于 0, 表 明 这 4 种 土地 
利用 类 型 的 面积 增加 ;林地 .中 覆盖 植被 和 未 利用 
土地 表现 为 单一 动态 度 均 小 于 0, 双 向 动态 度 均 大 
于 0, 表 明 这 3 种 土地 利用 类 型 的 面积 减少 。 结 合 六 
一 动态 度 和 双向 动态 度 发 现在 1990 一 2018 年 ,各 土 
地 利用 类 型 的 单一 动态 度 均 小 于 双向 动态 度 , 表 明 
在 此 期 间 绿 洲 内 部 各 土地 利用 类 型 的 相互 转化 较 
为 频繁 ,稳定 性 越 差 ,双向 动态 度 越 大 绿洲 土地 类 
型 越 容易 受到 外 界 因素 的 影响 。 
2.3 绿洲 植被 NDVI 变 化 特征 

归 一 化 植被 指数 作为 最 常用 的 植被 响应 指标 ， 
能 够 较为 精确 的 反映 地 表 不 同 覆 盖 程 度 植被 的 生 
长 态势 。 通 过 AreGIS 栅 格 计算 器 将 策 勒 绿洲 1990 
年 .2000 年 .2010 年 和 2018 年 的 4 期 策 勒 绿洲 NDVI 


表 2 1990—2018 年 策 勒 绿洲 土地 利用 动态 变化 率 


Tab. 2 Dynamic change rate of land use in Cele oasis from 1990 to 2018 


1990—2000 4 2000—2010 4E 2010—2018 年 1990—2018 年 
Ee | 变化 ”单一 动 ” 双 向 动 ” 变化 ”单一 动 ” 双 向 动 ” 变化 ”单一 动 ” 双 向 动 ” 变化 ”单一 动 ” 双向 动 
面积 km 态度 /% 态度 /% ”面积 km” 态度 /% 态度 /% ”面积 km” 态度 /% 态度 /% ”面积 km” 态度 /% 态度 /% 
农田 15.62 2.46 4.55 28.16 3.56 6.11 27.07 3.15 3.32 70.84 3.98 10.41 
林地 4.92 1.13 18.30 -41.55 -8.59 12.78 0.01 0.01 0.90 -36.62 -3.01 3.93 
高 覆盖 植被 22.18 3105 40.62 8.96 3.06 14.95 -10.50 -3.43 3.74 20.64 10.32 15.14 
pimi -40.00 -7.93 10.17 31.16 29.775 44.99 -8.53 -2.56 2.96 -17.37 -1.23 3.40 
低 覆 盖 植 被 7.14 1910 38.10 -2.93 -2.69 13.50 -L74 -2.73 3.57 2.47 2.36 9.10 
建设 用 地 2.69 5.00 16.19 -3.70 . -457 9.50 3.37 9.59 10.25 2.36 1.56 5.52 
未 利用 土地 -12.64 -1.52 409  -1975 -2.81 3.35 -9.73 | -240 2.68 -42.12 -1.81 2.06 
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求 均值 如 图 4a 所 示 , 由 图 可 看 出 ,多 年 来 该 绿洲 
NDVI 空间 分 布 表 现 为 中 心 大 于 四 周 ,NDVI 值 由 中 
心 向 四 周 递减 。 由 图 4 可 知 , 在 1990 一 2018 年 间 
策 勒 绿洲 NDVI 呈 现 出 不 同 程度 的 波动 ,波动 范围 
在 -0.46~0.60 之 间 ,红色 部 分 是 轻微 ,严重 退化 区 域 
占 比 为 27.56% ,其 中 空白 部 分 表示 NDVI 增 大 的 区 
域 , 占 整 个 绿洲 的 比重 为 72.44%。 结 合 策 勒 绿洲 土 
地 利用 类 型 空间 分 布 ,绿洲 植被 NDVI 增 加 的 土地 
利用 类 型 主要 为 未 利用 土地 ,绿洲 农田 内 部 NDVI 
退化 的 原因 可 以 概况 为 两 个 方面 :(1) NDVI 退 化 的 


量 ; 人 为 因素 包括 人 口 .粮食 产量 ,以 下 将 首先 对 各 
素 的 变化 特征 进行 分 析 ,然后 探讨 各 因素 对 绿洲 
土地 类 型 变化 的 影响 。 

2.4.1 驱动 因素 变化 特征 ”1960 一 2015 年 对 策 勒 绿 
洲 的 年 均 降 水 量 .年 均 温 .年 均 蒸 发 量 . 策 勒 河 径流 
量 、 人 口 .粮食 产量 采用 M-K 趋势 分 析 法 进行 研究 ， 
从 图 5 可 知 ,降水 量 在 1990 一 2000 年 间 处 于 显著 下 
降 趋势 ,2000 年 之 后 UF>0, 降 水 量 开 始 增加 ,但 增 
加 的 趋势 不 明显 。 年 均 温 从 1960 年 来 呈现 出 先 减 
少 后 增加 的 趋势 ,1997 年 之 后 UF 统计 量 均 超过 了 


区 域 与 绿洲 内 建设 用 地 的 扩 增 有 关系 ,是 由 于 在 农 
田 内 部 居民 建筑 用 地 的 扩 增 ,以 及 农田 灌溉 渠 系 的 
建立 和 向 四 周延 伸 , 随 着 建设 用 地 的 扩 增 ,绿洲 内 
SIUE SE NDVI REAR. (2) 1990 年 绿洲 内 部 分 布 林 
地 的 分 布 面积 较 大 ,到 2018 年 林地 面积 显著 缩减 ， 
1990 年 和 2018 年 绿洲 内 部 作物 和 植被 的 种 植 结 构 
不 一 样 ,而 据 前 人 研究 结果 “显示 :林地 的 NDVI 
值 大 于 农田 的 NDVI, 到 2018 年 绿洲 内 部 农田 区 域 
的 NDVI 值 呈现 出 降低 。(3) 1990 年 绿洲 内 部 农田 
种 植 主要 是 以 粮食 作物 为 主 ;2018 年 绿洲 内 部 种 植 
业主 要 以 经 济 作物 为 主 ,不 同 作 物 植 被 类 型 的 ND- 


VI 值 不 同 。 
2.4 土地 利用 的 驱动 因素 分 析 
策 勒 绿洲 是 典型 的 干旱 区 ,生态 系统 结构 简 


单 , 极 易 受 自然 和 人 类 活动 影响 ,对 气候 变化 响应 
敏感 ,小 幅度 的 气候 变化 都 将 对 该 绿洲 的 生态 环境 
产生 明显 的 影响 。 本 研究 主要 考虑 的 自然 因素 包 
括 : 策 勒 河 径 流量 、 年 均 降水 量 \ 年 均 温 、 年 均 蒸发 


显著 水 平 线 ,表明 温度 持续 增加 且 增 加 趋势 显著 。 
蒸发 量 从 1960 年 以 来 表现 出 先 减少 后 增加 的 趋势 ， 
在 1985 年 之 后 增加 趋势 线 超过 显著 性 水 平 ,说 明 获 
发 量 以 显著 的 趋势 持续 增加 。 因 此 ,说 明 策 勒 绿洲 
在 1990 年 之 后 气候 偏 暧 干 , 降 水 量 减少 ,蒸发 量 增 
加 ,绿洲 的 气候 变 和 干燥, 导致 地 表 水 的 大 量 莱 发 , 进 
而 导致 植被 .土壤 需 水 量 的 增加 ,最 终 影响 绿洲 植 
被 的 分 布 生长 状况 。 降 水 量 和 温度 的 变化 除了 直 
接 导 致 天 然 植被 面积 发 生变 化 ,还 通过 调节 土壤 中 
含水 量 间接 地 影响 绿洲 覆 被 的 空间 布局 。 而 人 口 
数量 在 1990 年 之 后 持续 增加 , 且 增加 趋势 显著 。 粮 
食 产 量 在 1990 年 以 后 增加 趋势 线 超过 置信 水 平 , 粮 
食 产量 明显 增加 。 人 口 增加 对 农田 用 地 的 需求 增 
加 ,农田 面积 发 生 扩 增 。 策 勒 河 径 流 1990 之 后 表现 
为 不 显著 的 增加 趋势 。 依 据 和 田 用 水 计划 ,从 策 勒 
河 提取 的 水 资源 约 有 82.1% 用 于 农业 灌溉 , 随 着 农 
田 面 积 和 人 口 的 持续 显著 增加 ,对策 勒 河 的 需求 量 
将 持续 增 大 。 现 有 水 资源 只 能 满足 饮水 `. 农 业 和 工 


N 
(b) 人 
NDVI 差 值 
mu -0.46~0.00 
C 0.00~0.60 0 3 km 


图 4 1990 一 2018 年 策 勒 绿洲 NDVI 平 均值 分 布 (a) 和 NDVI 差 值 变 化 (b) 


Fig.4 Change of mean NDVI distribution (a) and difference of NDVI (b) in Cele oasis from 1990 to 2018 
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图 5 各 驱动 因素 变化 趋势 


Fig. $ Variation trend of drivers 


业 的 需 水 , 当 农田 面积 增加 的 速度 超过 一 定 程度 
时 , 策 勒 河 径流 量 将 难以 满足 下 游 天 然 植 被 的 需 
水 ,进而 影响 草地 植被 的 分 布 与 数量 。 

2.4.2 驱动 因素 回归 分 析 为 排除 因子 之 间 的 相关 
性 问题 ,对 6 个 有 影响 因子 进行 主 成 分 分 析 , 并 依据 计 
算出 的 特征 值 和 累积 方差 贡献 率 的 大 小 确定 主 成 
分 数目 。 由 表 3 可 以 看 出 ,第 一 、 第 二 主 成 分 的 累计 
贡献 率 超过 90% ,说 明 前 2 个 主 成 分 已 经 能 充分 代 
表 原 始 驱 动因 子 来 反映 策 勒 绿洲 土地 利用 的 具体 
情况 。 进 一 步 计 算 6 个 驱动 因子 在 前 两 主 成 分 上 的 
载荷 ,发现 人 口 .粮食 产量 在 第 一 主 成 分 上 的 载荷 


较 大 ,说 明 第 一 主 成 分 代表 经 济 社会 因素 。 菊 发 
量 、 气 温 在 第 二 主 成 分 上 的 载 答 最 大 ,说 明 第 二 主 
成 分 仿 向 自然 因素 。 

再 运用 多 元 回归 分 析 法 ,计算 该 地 区 土地 利用 
类 型 变化 与 提取 出 的 2 个 主 成 分 之 间 的 函数 关系 。 
以 各 年 份 主 成 分 得 分 系数 作为 自 变 量 构 建 回归 算 
阵 , 进 而 建立 起 一 个 多 元 线性 回归 模型 ,在 计算 时 ， 
将 土地 利用 类 型 设 为 因 变量 Y; 主 成 分 一 为 自 变 量 
X1, 主 成 分 二 设 为 自 变 量 X2; 同 时 选用 检验 和 拟 
合 度 检验 确定 分 析 结 果 的 准确 性 ,回归 分 析 结 果 与 
回归 分 析 的 参数 估计 及 假设 检验 结果 见 表 4。 


表 3 土地 覆 被 变化 驱动 因子 主 成 分 分 析 


Tab.3 Principal component analysis of driving factors of land cover change 


因素 第 一 主 成 分 第 二 主 成 分 第 三 主 成 分 第 四 主 成 分 第 五 主 成 分 

径流 量 0.448 0.312 -0.203 -0.168 -0.752 
气温 0.290 0.692 0.344 0.512 0.221 
降水 0.440 0.043 0.411 -0.727 0.326 
蒸发 量 -0.325 0.598 -0.569 -0.369 0.246 
AH 0.451 -0.172 -0.479 0.190 0.468 
粮食 0.461 -0.190 -0.344 0.089 -0.004 
特征 值 4.327 1.140 0.336 0.194 0.003 
累计 贡献 率 /% 0.721 0.911 0.967 0.999 1 
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由 表 4 可 知 ,农田 用 地 变化 与 第 一 、 第 二 主 成 分 
均 表 现 为 正 相 关 性 , 且 显 著 性 水 平 较 高 ,其 中 与 代 
表 社 会 因素 的 第 一 主 成 分 的 相关 性 较 大 , 随 着 人 口 
的 增加 农田 用 地 明显 增加 。 建 设 用 地 与 第 一 ,第 二 
主 成 分 拟 合 出 的 模型 ,显著 水 平 最 高 (R=96.5% ) ， 
表现 为 与 第 一 主 成 分 呈正 相关 ,根据 主 成 分 代表 的 
含义 ,说 明 社 会 因素 驱动 建设 用 地 的 增加 , 随 着 社 
会 经 济 的 发 展 , 人 们 对 住址 基础 建设 等 需求 增加 ， 
建设 用 地 面积 随 之 扩 增 。 未 利用 土地 面积 变化 拟 
合 出 的 模型 显著 水 平 同 样 较 高 ,与 第 一 、 第 二 主 成 
分 均 为 负 相关 关系 ,首先 在 与 第 一 主 成 分 代表 的 社 
会 因素 方面 , 随 着 人 口 的 增加 ,绿洲 的 用 地 压力 变 
大 ,未 利用 土地 被 不 断 开 晨 ;在 与 第 二 主 成 分 代表 
的 自然 因素 方面 , 随 着 降水 量 和 温度 增加 ,促进 天 
然 植 被 的 生长 ,因而 未 利用 土地 面积 减少 。 其 余 土 
地 利用 类 型 与 两 主 成 分 拟 合 出 的 模型 显著 性 水 平 


积 不 断 减少 ,建设 用 地 面积 增加 。 绿 洲 各 土地 利用 
类 型 的 变化 状况 与 汉 瑞 英 等 关于 新 疆 南部 土地 
利用 变化 研究 结果 一 致 。 

(2) 策 勒 绿洲 土地 利用 类 型 发 生 明 显 转变 LT 
究 其 原因 ,主要 可 以 概况 为 人 为 自然、 政策 因素 这 
三 方面 。 在 人 为 因素 方面 :(1) 随 着 人 口 的 快速 增 
加 ,人 们 为 满足 日 常生 活 的 需求 ,过 度 的 采伐 树木 、 
放牧 破坏 草地 植被 和 开明 农田 ; (2) 策 勒 河上 游 建 
立 水 库 , 流 到 沙漠 中 的 水 量 会 进一步 减少 ,影响 绿 
洲 过 渡 带 以 胡杨 、 树 柳 、 骆 驼 刺 等 为 代表 的 自然 植 
被 乃至 整个 绿洲 的 存在 ;(3) 人 们 建立 渠 系 截流 河 
水 用 于 灌溉 农作物 而 间接 造成 自然 植被 的 大 面积 
衰退 或 死亡 的 ”。 在 自然 因素 方面 ,降水 、 径 流量 
起 的 作用 较为 明显 ,具体 表现 为 : 随 着 绿洲 面积 扩 
大 、 需 水 量 增 加 ,致使 流入 沙漠 的 水 量 减 少 ,下 游 植 
被 的 恢复 更 新 得 不 到 保证 , 且 下 游 水 量 的 减少 使 地 


较 低 ,说 明 这 些 土地 利用 类 型 面积 变化 与 驱动 因素 
的 相关 性 较 小 。 


表 4 回归 分 析 结 果 
Tab. 4 Regression analysis results 
土地 利用 类 型 回归 方程 显著 性 水 平 
农田 Y=0.429X1+0.282X2 R^-88.7496, F=11.82 


林地 Y--0.330X1-0.384X2 R^-63.9396, F-2.66 

高 覆盖 植被 Y-0.190X2-0.287X2 R^-25.0796, F-0.50 

中 覆盖 植被 y=-0.105X1+0.618X2  R=48.30%, F=1.40 

TIC si ELE Y-0.094X1-0.737X2 R=65.76%, F-2.88 

建设 用 地 Y-0.278X1-0.744X2 R^-96.5096, F241.40 

未 利用 土地 y=-0.426X1-0.161X2 — R'-81.4446, F=6.58 
3 讨论 


(1) 以 位 于 塔里木 盆地 南 缘 的 典型 脆弱 绿洲 为 
研究 区 ,基于 1990 年 .2000 年 .2010 年 .2018 年 策 勒 
绿洲 4 期 土地 利用 数据 和 同期 遥感 影像 数据 ,经 研 
究 发 现 该 绿洲 各 土地 利用 类 型 的 空间 分 布 位 置 TG 
积 均 发 生 显著 变化 ,不 同类 型 之 间 转 化 较为 频 发 ， 
转化 面积 占 到 总 面积 53.85% ,绿洲 动态 性 大 。 其 中 
农田 的 面积 显著 向 四 周 扩 增 ,尤其 向 东北 部 扩 增 最 
为 明显 。 林 地 的 空间 分 布 由 东北 -西北 -内 部 转移 ， 
且 面 积 不 断 减少 ,草地 的 分 布 区 域 不 断 向 绿洲 边缘 
缩减 ,再 者 根据 计算 的 林 草 植被 指数 同样 说 明 绿洲 
林 草 地 一 直 处 于 退化 的 状态 。 期 间 未 利用 土地 面 


下 水 补充 量 不 足 , 水 位 下 降 ,从 而 依靠 地 下 水 生长 
的 植被 衰退 .消亡 。 
(3) 绿洲 边缘 的 天 然 植被 并 非 都 是 耐 旱 植物 ， 
其 生命 的 维持 主要 是 靠 地 下 水 ,此 类 植被 一 旦 因 人 
类 破坏 造成 的 地 下 水 位 降低 而 衰亡 ,要 恢复 如 初 非 
常 难 ,最 终 影响 策 勒 绿洲 的 土地 利用 类 型 分 布 。 在 
政策 因素 方面 ,西部 大 开发 .塔里木 河 生 态 综合 治 
理工 程 .19 省 市 援 疆 等 政策 的 实施 ,其 中 天 津 对口 
文 援 策 勒 县 进行 了 沙漠 综 合 治理 、 建 立 灌区 试点 等 
计划 ,这 些 政策 促进 绿洲 经 济 不 断 发 展 ,人 口 数量 
持续 增加 ,为 满足 经 济 社会 发 展 需求 ,不 断 开 晨 土 
地 ,人 工 绿洲 面积 不 断 扩大 ,绿洲 植被 覆盖 发 生变 
化 ,直接 影响 绿洲 的 土地 利用 状况 。 综 合 而 言 , 通 
常 自 然 因 素 的 影响 在 短 时 间 内 造成 的 影响 是 很 小 
的 ,主要 限制 了 绿洲 的 低 覆 盖 度 植被 的 生存 ,其 主 
影响 是 长 期 累积 效应 。 而 人 为 因素 带 来 的 影响 
在 任何 时 期 内 都 是 非常 明显 ,人 文 因 素 是 绿洲 土地 
覆盖 变化 的 直接 驱动 力 , 对 土地 利用 变化 起 到 主导 
VERI", 


4 结论 


vd 


(1) 从 空间 分 布 上 来 看 , 策 勒 绿洲 土地 利用 类 
型 的 空间 分 布 发 生 明 显 的 转变 ,农田 分 布 面积 明显 
向 四 周 扩 增 ,其 中 向 东北 部 扩 增 最 为 显著 ,各 时 期 
林地 的 迁移 方向 为 东北 -西北 -内 部 ,草地 分 布 区 域 
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逐渐 向 绿洲 边缘 缩减 ,未 利用 土地 面积 不 断 减少 。 

(2) 从 时 间 尺 度 上 来 看 ,1990 一 2018 年 间 策 勒 
绿洲 共有 53.85% 的 面积 发 生 转 变 ,其 中 农田 、 居 民 
建设 用 地 面积 增加 ,林地 、 未 利用 土地 .高 覆盖 草 
地 .中 覆盖 草地 和 低 履 盖 草 地 面积 减少 。 各 土地 利 
用 类 型 的 动态 度 变 化 表现 为 :农田 高 覆盖 植被 AR 
窗 盖 植被 和 建设 用 地 的 单一 动态 度 和 双向 动态 度 
均 大 于 0, 表 明 这 4 种 土地 利用 类 型 的 面积 增加 ; 林 
地 .中 和 覆盖 植被 和 未 利用 土地 表现 为 单一 动态 度 均 
小 于 0, 双 向 动态 度 均 大 于 0, 林 草 植被 指数 小 于 0。 

(3) 从 NDVI 的 分 布 上 来 看 ,表现 为 中 心 大 于 四 
周 ,NDVI 值 由 中 心 向 四 周 递减 。 从 NDVI 的 变化 方 
面 看 ,表现 为 NDVI 增 大 的 区 域 (72.44% ) 主 要 集中 
在 绿洲 边缘 沙漠 分 布 区 , NDVI 减少 的 区 域 
(27.56%) 主要 集中 在 绿洲 内 部 农田 分 布 区 。 

(4) 各 驱动 力 因 素 对 土地 利用 类 型 的 影响 结果 
为 :建设 用 地 农田 和 未 利用 土地 受 各 因素 的 影响 
最 显著 ,农田 和 建设 用 地 面积 变化 与 人 为 因素 呈正 
相关 性 ,而 未 利用 土地 面积 变化 与 人 为 因素 和 自然 
因素 均 表 现 为 负 相 关 性 。 

(5) 干旱 区 绿洲 生态 系统 结构 简单 , 极 易 受 自 
然 和 人 类 变化 的 影响 , 且 反 应 显著 。 本 文 从 土地 利 
用 转移 矩阵 .土地 利用 动态 度 . 林 草 植被 指数 3 个 指 
标定 量 定性 的 探讨 了 策 勒 绿洲 的 土地 利用 /土地 覆 
被 变化 ,并 对 导致 其 变化 的 驱动 因素 进行 了 分 析 ， 
为 合理 开发 利用 绿洲 提供 科学 依据 。 
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Abstract: The characteristics of land use/cover change of the Cele oasis, in the southern margin of Tarim Basin, 


from 1990 to 2018 were explored and attributed based on four phases of land use and remote sensing images, as 


well as meteorological, hydrological, and socio-economic data. We found that the land use types in the Cele oasis 


changed substantially: The farmland area expanded to the surrounding area, the forest land distribution area 
shifted substantially, and the coverage of grassland changed frequently. Overall, 138.41 km" (53.8594) of the total 


area was transformed. In descending order, the single dynamic attitude was high covered grassland>farmland> 


low cover grassland»construction land?medium covered»grassland unused land»forest land, and in descending 


order, the double directional trend was high covered grassland»farmland-low cover grassland»construction land 


forest land>medium covered»grassland unused land. The results also showed that the vegetation index of the 


green continent was less than 0, which indicates that forest land and grassland have been in a state of degradation. 


The driving forces of land use change in the oasis were population, precipitation, and runoff, based on the factors 


influencing surface cover in the oasis. 


Keywords: desert oasis; Cele oasis; land coverage; spatio-temporal variability; NDVI; driving force 


